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ABSTRACT 

This  is  a  report  on  the  boundary  layer  aspects  of  the  NPS 
participation  in  CEWCOM-78.   The  primary  purpose  of  the  experi- 
ment was  to  determine  how  representative  San  Nicolas  Island  is 
of  an  open  ocean  marine  boundary  layer  and  to  examine  the  vali- 
dity of  boundary  layer  measurements  at  the  NRL  tower  on  the  MW 
tip  of  the  island.   Under  favorable  wind  conditions  (NW)  the 
turbulence  and  profile  structure  of  the  boundary  layer  near  SNI 
was  characteristic  of  typical  marine  conditions.   A  comparison 
of  simultaneous  measurements  at  the  NRL  tower  and  the  R/V 
Acania  indicated  considerable  shoreline  influence  on  the  velo- 
city fluctuations  (U*  or  e)    and  the  mean  wind  speed  (U)  but 

2 

essentially  no  influence  on  temperature  fluctuations  (C„  ) . 

Using  the  bulk  method  to  calculate  T^  and  £  from  the  Acania 

2 

data,  the  actual  measurements  of  C  "  could  be  predicted  to 

within  about  a  factor  of  two. 
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I .  SUMMARY 

A.  Introduction 

This  is  a  report  on  atmospheric  measurements  made  by 
the  Naval  Postgraduate  School  Environmental  Physics  Group  aboard 
the  R/V  Acania  in  the  vicinity  of  San  Nicolas  Island  (SNI)  in 
May  of  1978.   The  primary  goal  of  these  measurements  was  to  exa- 
mine the  open  ocean  "representativeness"  of  SNI  and  to  evaluate 
the  validity  of  measurements  made  at  the  NRL  tower  site  on  the 
north  west  tip  of  the  island.   This  report  will  focus  on  the 
turbulence  and  boundary  layer  data,  leaving  the  aerosol  evalua- 
tion for  a  later  report.   In  addition,  the  NPS  group  provided 
direct  micrometeorological  support  for  optic  experiments  and  a 
rather  conclusive  study  of  the  bulk  method  scaling  law  predic- 
tions of  turbulence  parameters  in  the  surface  layer. 

B.  Conclusions 

How  representative  SNI  is  of  the  marine  condition  is 
more  of  an  aerosol  question  than  a  turbulence  question.   How- 
ever, the  turbulence  aspect  is  important.   During  CEWCOM-76  it 

was  found  that  coastal  areas  exhibited  diurnal  variations  of 

2 

temperature  structure  function,  CT  ,  characteristic  of  overland 

2 
sites  (minima  in  C„   at  sunrise  and  sunset)  whereas  open  ocean 

areas  exhibited  almost  no  diurnal  variation.   Under  the  W-NW 

wind  conditions  that  predominated  during  the  turbulence  evalua- 

'T 


2 

tion  periods  of  CEWCOM-78,  the  C„   measurements  near  SNI  showed 


no  obvious  diurnal  variation. 

The  Naval  Research  Laboratory  (NRL)  tower  site  mea- 
surement made  by  NPS  personnel  using  identical  equipment  to  that 


being  employed  on  the  R/V  Acania  have  been  compared  to  simul- 
taneous shipboard  measurements.   For  data  taken  when  the  Acania 

was  in  the  immediate  vicinity  of  the  tower  site  (primarily 

2 
anchored  within  .3  miles)  the  C_   comparison  showed  excellent 

agreement.   Neglecting  a  few  low  wind  speed  cases,  the  average 
disagreement  was  only  7%  for  23  periods  with  a  single  measure- 
ment standard  deviation  of  64%.   Given  the  combined  measurement 
uncertainty  of  about  30%  and  the  uncertainty  introduced  by  the 
stochastic  nature  of  atmospheric  turbulence,  it  may  be  quite 

difficult  to  do  significantly  better.   Although  only  a  few 

2 

periods  were  available,  C„   measurements  made  at  the  NRL  tower 

compared  fairly  well  with  shipboard  measurements  made  when  the 
Acania  was  30  to  50  miles  upwind  of  the  island.   This  indicates 
that  SNI  is  in  a  region  of  good  horizontal  homogeneity  under 
NW  wind  conditions. 

Comparisons  of  wind  speed,  U,  and  the  rate  of  disipa- 

tion  of  turbulence  kinetic  energy,  e  ,  were  not  nearly  as  favor- 

2 

able  as  the  C™'  comparison.   The  values  of  e  were  used  to  cal- 
culate the  friction  velocity  U*  (the  surface  stress  is  propor- 

2 
tional  to  U*  ) .   On  average,  the  tower  measurements  of  UA  were 

2.5  times  greater  than  the  ship  measurements  with  a  standard  de- 
viation of  93%.   The  tower  measurements  of  wind  speed  (at  Z=ll 
meters)  were,  on  average,  16%  lower  than  the  ship  measurements 
with  a  standard  deviation  of  10%.   The  lower  wind  speed  and 
higher  surface  stress  at  the  tower  is  a  result  of  the  increased 
drag  imposed  by  the  surf  and  land.   This  means  that  neither 


turbulence  nor  profile  measurements  at  the  tower  can  be  used  to 
determine  the  atmospheric  stability  and  Monin-Obukhov  scaling 

parameters  over  the  immediate  ocean  area. 

2 
The  estimation  of  Cm   (as  well  as  e)  using  Monin- 

Obukhov  scaling  parameters  not  only  requires  a  validation  of 

2 

the  CT   parameterization  formulae,  but  also  requires  a  prac- 
tical method  of  obtaining  the  scaling  parameters.   Employing 
only  four  physical  quantities  (sea  surface  temperature  and 
wind  speed,  air  temperature  and  relative  humidity  at  some  re- 
ference height  above  the  sea  surface)  the  bulk  method  is  not 
only  the  simplist  but  is  also  the  least  demanding  in  terms  of 

accuracy.   Using  data  from  several  cruises,  the  NPS  group  has 

2 
shown  that  Wyngaard  et  al.'s  (1971)  Cm   parameterization  is  valid 

over  the  ocean  and  that  the  bulk  method  provides  an  excellent 
rendition  of  the  scaling  parameters  [Davidson  et  al.  (1978)]. 
Included  in  this  report  is  a  comparison  of  bulk  predictions  with 
observations  obtained  during  CEWCOM-78,  demonstrating  the  appli- 
cability of  this  technique. 

C.    General  Comments  and  Recommendations 

1.    Obviously  the  NFL  tower  site  location  is  most 
suitable  for  NW  wind  conditions.   During  CEWCOM-78  we  did  find 
good  tower  data  for  winds  from  240  to  320  degrees  (this  does 
not  necessarily  exclude  other  wind  directions).   However,  there 
is  evidence  of  land  influence  under  light  wind  conditions.   In 
view  of  this,  it  would  be  prudent  to  limit  the  tower  data  to 
wind  speeds  greater  than  2.5  m/sec. 


2.  We  feel  that  atmospheric  stability  and  Monin- 
Obukhov  scaling  parameters  should  be  calculated  using  the  bulk 
method  with  the  updated  coefficients  and  techniques  given  in 

the  text.   This  will  require  supplementing  the  tower  measurements 
with  a  sea  surface  or  bulk  ocean  temperature  measurement.   Also, 
the  tower  measured  wind  speed  should  be  increased  to  account  for 
the  surface  drag  effects. 

3.  Given  the  importance  of  the  humidity  contribution 

2 

to  CN  ,  the  temperature-humidity  covariance,  C„  ,  should  be 

measured  at  the  tower  site.   Based  on  bulk  estimates,  during 

2        2 
CEWCOM-78  the  average  relative  contribution  of  CT   and  C„   was 

2 
70%,  the  C    contribution  was  24%  and  the  C    contribution  was 

6%. 

II.   SHIP  MOVEMENTS 

The  primary  movements  of  the  R/V  Acania  are  shown  in  Fi- 
gures 1  and  2.   Anchorage  locations  at  SNI  are  shown  in  Figure 
3.   A  summary  of  data  periods  relevant  to  SNI  evaluations  is 
given  in  Table  I.   Periods  of  running  downwind  and  periods  in- 
side the  Channel  Islands  have  been  excluded. 
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III.  THEORETICAL  SUMMARY 
A.    Definitions 

The  optically  relevant  turbulence  quantities  measured 

2 
were  temperature  structure  function  parameter,  CT  ,  and  the  rate 

of  dissipation  of  turbulent  kinetic  energy,  z.      The  refractive 

2  2 

index  structure  function  parameter, CN  ,  is  related  to  C_   by 

CN2  =  (79  x  10"6  P/T) (CT2  +  .11  CTq  +  .0032  Cq2)     (1) 

where  C_   is  the  temperature  -  water  vapor  covariance  parameter 

2 

and  C    is  the  water  vapor  structure  function  parameter.   The 

microscale  of  turbulent  functions,  n ,  is  related  to  £  by 

n  =  (v3/e)1/4  (2) 

where  v  is  the  kinematic  viscosity  of  air. 
R.    Monin-Obukhov  Similarity  (MOS) 

In  the  atmospheric  surface  layer  (Z  =  50  meters)  C„ 
and  £  can  be  calculated  from  scaling  parameters  that  are  related 
to  more  easily  measured  atmospheric  properties  [see  Haugen  (19731 
for  a  review] .   The  appropriate  MOS  formulae  are 


2  _  m  2  „-2/3 
*T 


C„      =  T*  Z~"/0    f  (£)  (3) 


£  =  (U*3/Kz)E(£)  (4) 

Where  T^  and  U*  are  the  potential  temperature  and  wind  speed 
scaling  parameters,  f(£)  and  E(£)  are  dimensionless  MOS  func- 
tions (Appendix  A) ,  Z  is  the  vertical  coordinate  and  K  is  Von 
Karmon's  constant  (K  =  .35).   The  MOS  dimensionless  stability 
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parameter,  £,  is  related  to  the  Monin-Obukhov  length,  L, 


KaZ 


(T*  +  .61  Tq*) 


K   =  Z/L  =  ^  , (5) 


U  2 


'* 

where  q*  is  the  water  vapor  mixing  ratio  scaling  parameter,  T 
is  the  temperature  and  g  is  the  acceleration  of  gravity. 

The  vertical  profile  of  the  mean  quantity  X  (where 
X  =  T,  q,  U)  can  be  represented  by 

x* 
X(Z)  =  X(0)  +  ^-^    (In  Z/Zqx  -  *x<5>)  (6) 

where  Z    is  the  roughness  length  for  X  and  ^V(C)  is  the  profile 

function  (Appendix  A) .   The  value  of  K  is  chosen  so  that  a  =1. 

We  have  assumed  am  =  a   =1.35. 

T     q 

C.    Bulk  Method 

2 
One  cannot  calculate  CT   and  £  from  Equations  3  and  4 

until  one  first  obtains  values  for  T^ ,  U*  and  £.   The  bulk  method 

of  determining  the  MOS  scaling  parameters  is  based  upon  relating 

X*  to  the  air-sea  X  difference  (  A  X)  through  the  drag  coefficient 

cx- 


,X1/2(X(Z)  -  X(0))  =  cxly 


X*  =  cv  '     (X(Z)  -  X(0))  =  cv  '      AX  (7) 


Using  Equation  6  we  can  define  the  neutral  value  (£  =  0)  of  the 
drag  coefficient,  c  N 

c   1/2  =     axK 

XN       In  Z/Z  (b) 

ox 

and  relate  cx„  to  cy 

^XN 

°X     ,-,     ,  -*-l         1/2       "> 


(1  -  (axK)   cXN    *  (£)) 
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Using  Equation  5  and  Equation  7  we  can  calculate 

(i  -  K^c^^m)2 

(1  -  (aTK)  1cTN-L/^T(C)) 
with 

KaZ       CTN1/2(AT    +     •18Ac3) 

r    =  K<?z  _i£! (ID 

°    T  c   U2 

°UNU 

where  A  T  is  the  air-sea  potential  temperature  difference  (°C) , 
Aq  is  the  air-sea  water  vapor  mixing  ratio  difference  (gm/kg)  and 

U  is  the  wind  speed  (m/sec) .   c,.«  varies  with  wind  speed  but  is 

-3 
well  approximated  by  cUN  =  1.3  x  10   .   Based  upon  Davidson  et 

_  3 
al.'s  (1978)  work,  a  good  value  for  c _„  is  cTN  =  1.3  x  10 

The  actual  bulk  method  process  goes  as  follows: 

1)  From  U  calculate  c„N 

2)  From  U,  AT,  Aq,  cUN  calculate  £   (Equation  11) 

3)  From  £  ,  solve  Equation  10  to  find  £ 

4)  From  £  and  cTN  calculate  cT  (Equation  9) 

5)  From  AT  and  cT  calculate  T+  (Equation  7) 

2 

6)  From  T*  and  £  calculate  CT   (Equation  3) 

The  process  can  be  greatly  simplified  by  ignoring  the  wind  speed 
dependence  of  cn„.   In  this  case  (for  z  =  10  meters  and  T=15°C) 

£Q  =  3.3  (AT  +  .18Aq)/U2 

We  have  solved  Equation  10  for  this  case,  allowing  a  simple  alge- 
braic relation  between  £   and  £. 
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e  -  50  a  -  -o3?0'4)  ?o<0  (13a) 

5  =  ?o(l  +  .18S0'8  +  .13?03)       ?o>0       (13b) 

This  leads  to  the  simplified  bulk  method  process: 

1)  From  U,  AT  and  Aq  calculate  £   (Equation  12) 

2)  From  £   calculate  £  (Equation  13) 

3)  From  £  and  c_N  calculate  cT  (Equation  9) 

4)  From  AT  and  cT  calculate  T*  (Equation  7) 

2 

5)  From  T*  and  £  calculate  C™   (Equation  3) 

2 

D.    Application  of  Bulk  Method  to  CN 

2 
In  order  to  calculate  CXT   from  Equation  1,  one  must 

2  2 

have  available  estimations  of  C^  ,  C    and  C   .   Since  the  bulk 

method  calculations  have  given  us  T* ,  q*  and  £,  let  us  suppose 

2 

that  C    and  CT   can  be  calculated  from  the  MOS  analagous  form 

for  Equation  3 

Cq2  =  q*2  Z"2/3  f  (£)  (14a) 

CTq  =  T^q+Z"2/3  f(£)  (14b) 

From  Equation  1,  we  now  have 

CN2  =  (79  x  10"6P/T2)2(T*2  +  .11  T*q*  +  . 0032q*2 ) Z~2/3f ( £ )   (15) 

2 

In  Equations  14a  and  14b  we  have  assumed  that  C    and  C™   obey 

Monin-Obukhov  similarity  and  they  have  the  same  dimensionless 

2 

structure  function  parameter  (f(£))  as  CT  . 

IV.   DATA 

A.    Surface  Layer  Data 

The  shipboard  turbulence,  mean  and  MOS  scaling  para- 
meter data  are  shown  in  Figures  4  through  12.   The  turbulence 
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Figure  4.  A/V  Acania  Atmospheric 
Surface  Layer  Data  for 
5/14/78. 
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Figure  5.  R/V  Acania  Atmospheric 
Surface  Layer  Data  for 
5/15/78. 
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Figure  6.  R/v  Acania  Atmospheric 
Surface  Layer  Data  for 
5/18/78. 
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Figure  7.  R/V  Acania  Atmospheric 
Surface  Layer  Data  for 
5/19/78. 
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Figure  8.  R/V  Acania  Atmospheric 
Surface  Layer  Data  for 
5/20/78. 
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Figure  9.  R/V  Acania  Atmospheric 
Surface  Layer  Data  for 
5/21/78. 
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Figure  10.  R/V  Acania  Atmospheric 
Surface  Layer  Data  for 
5/22/78. 
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Figure  11.  R/V  Acania  Atmospheric 
Surface  Layer  Data  for 
5/23/78. 
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figures  include  both  the  measured  values  and  the  values  predicted 
on  the  basis  of  the  bulk  method  MOS  parameters.   Data  tables  are 
given  in  Appendix  B. 

B.    NRL  Tower  Site  Evaluation 

Although  there  is  evidence  that  shore  based  measure- 

2 
ments  of  C_   can  be  poorly  correlated  with  measurements  made 

over  the  ocean  [Davidson  et  al.  (1976)],  there  is  still  hope 

that  the  NRL  tower  can  provide  meaningful  data.   The  primary 

source  of  this  optimism  is  the  excellent  location  of  the  tower 

on  a  well  exposed  point  of  the  island.   In  order  to  compare  NRL 

tower  measurements  with  R/V  Acania  measurements,  NPS  personnel 

2 
installed  and  operated  standard  NPS  C™   and  z    equipment  on  the 

NRL  tower.   In  this  way,  measurements  with  identical  equipment 

and  procedures  could  be  compared,  thus  eliminating  (or  reducing) 

one  source  of  uncertainty.   We  did  take  the  liberty  of  using  the 

wind  speeds  measured  by  NRL  equipment. 

A  table  of  the  NPS  measurements  from  the  NRL  tower  is 

given  in  Appendix  C.   Table  II  is  a  compilation  of  data  for  those 

time  periods  when  simultaneous  SNI  and  Acania  measurements  are 

available.   The  data  is  further  restricted  to  time  periods  when 

the  Acania  was  either  at  anchor  or  underway  immediately  upwind 

of  the  tower  site  on  the  NW  tip  of  SNI.   Included  at  the  end  of 

the  table  is  a  comparison  when  the  Acania  was  30  miles  to  50 

2 
miles  upwind  of  SNI  (5/20-1515  to  5/21-1045).   The  CT   and  U 

data  from  the  Acania  are  Z  =  10  meter  equivalent  values.   The 

2 
SNI  CT   and  U  data  are  from  the  Z  =  11.4  meter  level.   Since 

the  E  data  is  subject  to  greater  statistical  scatter,  the 
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TABLE  II.  Comparison  of  turbulence  (Ct^  and 
(U)  measurements  taken  simultaneously  at  the 
the  R/V  Acania. 


U*)  and  mean  wind  speed 
NRL  tower  and  aboard 


DATE 

TIME 

<j)  ,DEG 

cT2,io  3 

°C  /M"7 

f 

U* 

M/SEC 

U 

M/SEC 

SNI 

ACANIA 

SNI  • 

ACANIA 

SNI 

ACANIA 

5/14 

1315 

254 

5.76 

.40 

.19 

.096 

2.0 

2.1 

1515 

240 

4.40 

.25 

.12 

.094 

1.5 

2.0 

1745 

2  57 

.10 

.49 

.17 

.077 

3.3 

2.4 

1845 

241 

4.63 

.19 

.48 

.080 

2.4 

2.8 

1915 

240 

.51 

.26 

.55 

.13 

3.7 

4.4 

2015 

244 

.17 

.  29 

.50 

.14 

4.0 

4.5 

5/15 

0845 

276 

.33 

.38 

1.0 

.30 

9.0 

10.1 

0900 

27  5 

.22 

.27 

1.0 

.34 

10.5 

10.7 

0945 

279 

.45 

.32 

.86 

.33 

10.0 

11.6 

1000 

277 

.32 

.37 

1.0 

.32 

9.5 

11.4 

1515 

321 

1.49 

1.0  5 

.71 

.35 

11.7 

13.4 

1615 

276 

.148 

1.17 

.64 

.41 

11.9 

14.7 

1915 

279 

1.13 

1.44 

.73 

.45 

14.4 

17.3 

2000 

271 

1.72 

1.6  5 

.80 

.50 

14.3 

19.6 

5/18 

2045 

308 

.96 

1.01 

.19 

.10 

2.4 

3.0 

2200 

307 

.57 

2.87 

.13 

.15 

2.6 

6.3 

5/19 

1045 

300 

.392 

1.23 

.33 

.14 

3.4 

3.8 

1145 

290 

3.19 

4.17 

.43 

.15 

3.9 

4.9 

1315 

287 

5.87 

4.97 

.45 

.16 

4.1 

5.3 

1345 

285 

8.15 

5.26 

.58 

.17 

4.1 

5.3 

1845 

269 

5.99 

7.82 

.46 

.24 

5.7 

7.2 

2115 

217 

4.33 

6.63 

.54 

.23 

5.5 

7.1 

5/20 

0945 

326 

5.11 

8.00 

.60 

.25 

5.7 

7.6 

1045 

304 

4.74 

7.24 

.60 

.30 

5.7 

8.4 

1115 

297 

7.26 

4.25 

.76 

.30 

5.4 

7.6 

5/20 

1515 

274 

6.62 

5.34 

.80 

.35 

6.7 

8.6 

1*45 

302 

3.20 

5.16 

.90 

.42 

7.1 

10.6 

1945 

296 

3.67 

3.17 

.80 

.45 

7.6 

10.6 

2045 

30  5 

5.18 

1.77 

.76 

.43 

7.3 

10.3 

5/21 

1015 

398 

5.8  3 

2.^7 

.73 

.17 

5.8 

7.2 

1045 

304 

5.53 

1.90 

.60 

.16 
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r 

!8 

multi-level  measurements  were  combined  (using  Equation  4)  into 

2 
a  single  U*  value.   Only  the  Z  =  11.4  meter  SNI  C™   data  was 

used  for  the  comparison  since  those  measurements  were  available 

for  the  entire  period  from  5/14  to  5/21. 

2 

The  CT   comparison  (Figure  12)  indicates  very  reason- 
able agreement  between  the  tower  and  shipboard  measurements.   A 
few  cases  of  greater  disagreement  occurred  during  light  winds 
(U~2  m/sec) .   The  average  ratio  of  the  SNI  to  Acania  values 
(Table  III),  excluding  the  light  wind  cases,  indicates  only  a 
7%  average  disagreement  with  a  standard  deviation  of  64%.   Given 
a  single  measurement  error  of  about  20%,  the  combined  instrumen- 
tal error  for  this  comparison  is  about  30%.   On  the  last  two 

days  of  the  direct  comparison  (5/19  and  5/20)  we  also  had  avail- 

2 

able  CT   measurements  at  a  second  level  on  the  NRL  tower  Z=17.5 

meters) .   When  converted  to  Z=10  meter  equivalent,  the  level  2 

values  did  not  agree  as  well  with  the  Acania  values.   For  this 

2 
period,  a  level  1  and  2  combined  C™   would  have  a  20%  average  dis- 
agreement with  the  Acania  values,  still  a  respectable  result. 

The  velocity  data  shows  significantly  greater  shore- 
line influence  (Table  III) .   The  velocity  scaling  parameter 
comparison  indicates  U*  values  2.5  times  greater  at  the  tower 
than  at  the  Acania  (Figure  12) ,  implying  a  considerably  greater 
surface  drag  immediately  upwind  of  the  tower.   This  increased 
drag  leads  to  16%  lower  wind  speeds  at  the  tower  (Figure  12)  for 

the  Z  =  11.4  meter  level.   The  average  velocity  drag  coefficient 

-2 
measured  at  the  tower  was  crTlvT  =  1.0  x  10   ,  considerably  greater 

U  IN 

than  the  typical  ocean  value  of  1.3  x  10 
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Figure    12. 


Comparison  of  Turbulence  (CT   and  U*)  and  Mean 
Wind  (U)  Data  from  the  R/V  Acania  with  NPS  Mea- 
surements at  NRL  Tower  on  San  Nicolas  Island. 
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TABLE  III. 


RATIO  OF  SNI  TO  ACANIA  MEASUREMENTS 
(r)  FOR  CT2,  U*  AND  U 


CT 

u* 

u 

<r> 

1.07 

2.48 

.84 

a 

.64 

.93 

.10 

N 

23.0 

25.0 

25.0 

a\TN* 

.13 

.19 

.02 

We  also  compared  SNI  CT   measurements  to  Acania  values 

when  the  ship  was  30  miles  to  50  miles  upwind  of  the  island  (Fi- 

2 
gure  13)  .   Since  only  a  few  C"T   values  were  available  from  the 

SNI  measurements;  the  comparison  is  rather  incomplete  but  it  is 

encouraging  that  no  large  disagreements  were  found. 
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FIGURE  13.   Comparison  of  CT   Data  from  the  R/V  Acania 
with  NPS  Measurements  at  the  NRL  Tower  on 
San  Nicolas  Island.  For  this  data,  the  Acania 
was  30  miles  to  50  miles  upwind  of  SNI. 
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C.    Acoustic  Sounder  Data 

The  NPS  acoustic  echosounder  (Aeroenvironment  Model 
300)  was  mounted  aboard  the  Acania  and  provided  inversion  height 
and  plume  structure  data  throughout  the  cruise.   A  tabulation  of 
the  sounder  data  is  given  in  Appendix  D. 
V.    CONCLUSIONS 

A.    Bulk  Method 

A  brief  examination  of  Figures  4  through  11  will  re- 
veal that  the  bulk  method  is  a  fairly  good  predictor  of  the 

2 

measured  values  of  C-   and  e .   A  log  average  of  the  ratio  of 

2 
predicted  CT   to  measured  yields  a  single  measurement  standard 

deviation  of  120%   for  CEWCOM-78.   In  other  words,  the  bulk 

method  predicted  the  measured  value  to  within  roughly  a  factor 

of  2  for  a  single  half-hour  averaging  period.   A  considerable 

portion  of  this  scatter  is  due  to  measurement  problems  (such 

2 
as  salt  contamination  of  C    sensors)  discussed  in  Davidson 

et  al.  (1978).   Figure  14  [from  Davidson  et  al.  (1978)]  indi- 
cates that  excellent  average  agreement  between  the  bulk  method 

2 

and  measured  values  of  CT  . 

Included  in  the  boundary  layer  data  (Appendix  B)  is 

2 
a  calculation  of  CN   (Equation  15)  from  bulk  values  of  T* ,  q* , 

and  E,.      The  average  relative  magnitude  of  the  terms  contribut- 
ing to  CN2  were:  C^2    (70%),  CTq  (24%)  and  C  2  (6%).   This 
agrees  well  with  the  direct  measurements  during  the  BOMEX  ex- 
periments [Friehe  et  al.  (1975)]. 
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FIGURE    15. 


Dimensionless  temperature  structure  function 
(CT2  z2/3/T^2  =  f(£))  vs.  atmospheric  stability 
(£) .   The  circled  points  are  data  and  the  solid 
curve  is  Wyngaard  et  al's  (1971)  formulae  for  f(£) 
The  bulk  method  was  used  to  calculate  T*  and  £. 


B.    SNI  Boundary  Layer  Measurements  at  the  NRL  Tower 

Comparison  of  the  NPS  measurements  on  the  Acania  and 
the  NRL  tower  (Figure  12)  used  to  evaluate  the  tower  site  data 
were  summarized  in  Tables  II  and  III.   During  the  periods  of 
the  comparison  the  wind  direction  was  most  favorable  (westerly 
to  northerly)  so  be  advised  that  our  conclusions  should  be  con- 
fined to  those  limits.   Based  on  this  data  and  the  discussion  of 
Section  IV,  we  conclude  the  following  about  the  NRL  tower  site: 
1)   The  data  is  subject  to  land  influence  effects  for 
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wind  speeds  less  than  about  2.5  m/sec. 

*T 


2 
2)   The  C™   data  is  quite  representative  of  the  marine 


boundary  layer. 

3)  The  U*  (or  e)  data  is  not  representative  of  the 
marine  boundary  layer,  presumably  due  to  increased  surface  drag 
upwind  of  the  tower. 

4)  The  mean  wind  speed  (U)  is  about  15%  less  than  the 
marine  boundary  layer  value  at  Z  =  11.4  meters. 

5)  The  marine  boundary  layer  value  of  U*  cannot  be 
obtained  from  wind  speed  profiles  at  the  tower.  This  follows 
from  3  and  4. 

In  addition  to  these  conclusions,  we  offer  the  spe- 
culation that  T*  and  q^  cannot  be  correctly  inferred  from  tower 
measurements  of  temperature  and  humidity  profiles .   This  opin- 
ion is  based  upon  the  known  strong  interaction  between  velo- 
city structure  and  scalar  structure.   On  the  other  hand,  since 

2 

the  C™   data  seems  to  be  unaffected  at  the  tower,  it  is  quite 

possible  that  the  temperature  and  humidity  profiles  are  simi- 
larly unaffected.   Perhaps  this  question  can  be  resolved  by  a 
comparison  of  NRL ' s  profile  data  and  the  NPS  R/V  Acania  data. 
Alternatively,  one  could  use  the  bulk  method  with  the  tower  data 
(requiring  the  addition  of  a  suitable  sea  surface  sensor)  as  has 
already  been  suggested  by  Carl  Friehe. 

C.    SNI  as  a  Representative  Boundary  Layer 

From  a  turbulence  and  boundary  layer  point  of  view, 
SNI  seemed  to  be  a  good  example  of  a  typical  open  ocean  marine 
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boundary  layer  during  CEWCOM-78  (5/14-5/25) .   This  is  based 
primarily  on  the  lack  of  diurnal  variation  of  boundary  layer 

parameters  which  is  typical  of  open  ocean  conditions.   A  com- 

2 
parison  of  C-   values  (admittedly  limited)  showed  a  reasonble 

correlation  between  SNI  values  and  values  measured  30  miles  to 
50  miles  upwind.   Obviously  this  conclusion  is  limited  to  fa- 
vorable wind  directions. 
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APPENDIX  A 
MOS  STABILITY  FUNCTIONS 


The  forms  of  the  mean  gradient  functions  (Businger  et  al.,  1971) 

(JjyU)  -  (1  -  155)"1/4  K<0 

4>D(5)  =  (1  +  4.7?)  C>o 

(J)T(C)  =  (1  -  19C)"1/3  ?<0 

4)T(?)  =  (1  +  6.4?)  C>0 

The  mean  profile  function 

2 

<M?)    =    l-MO-31n<MO+21n(^ll£i)    +    2tan~1(J)  (  O  -ir/2  +  ln  ( 1+<t>n  (  g)  ) 


The  dimensionless  velocity  dissipation  function  (Wyngaard  and 
Cote,  1971) 

o/3  3/9 

Eu(5)  =  (1  +  ,51V/    )    '*-  ?<0 

EU(C)  =  (1  +  2.5?2/3)3/2  C>0 

The  dimensionless  temperature  structure  function  parameter 
(Wyngaard  et  al.,  1971). 

fT(C)  =  4.9(1  -  7C)'2/3  C<0 

fT(C)  =  4.9(1  +  2.4£2/3)  £>0 
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APPENDIX  B 


Marine  boundary  layer  evaluation  from  R/V  Acania  measure- 
ments during  CEWCOM-7  8  -  bulk  data,  MOS  scaling  parameters,  tur- 

2 

bulence  data  and  bulk  calculation  of  C„  .   The  bulk  and  turbu- 


lence data  are  Z  =  10  meter  equivalent  values.   The  MOS  scaling 

N 


2 

parameters  are  calculated  by  the  bulk  method.   The  CM   values 


are  bulk  estimates  from  Equation  15. 
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APPENDIX  C 

NPS  data  from  the  tower  site  at  SNI .   Level  1  is  a  Z=11.4 
meters  and  level  2  is  at  Z=17.5  meters.   The  wind  direction 
(<J>)  is  from  the  Acania  measurements.   The  drag  coefficient 
(c,-.)  is  from  the  NRL  tower  U*  data. 
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APPENDIX  D 
ACOUSTIC  ECHOSOUNDER  RESULTS  FROM  R/V  ACANIA 


P:  Surface  Plume  Maximum  Height  (m) 

Inl:  Lowest  Inversion  Height 

In2:  Second  Lowest  Inversion  (m) 

In3:  Third  Lowest  Inversion  (m) 

In4:  Fourth  Lowest  Inversion  (m) 

In5:  Fifth  Lowest  Inversion  (m) 

In6:  Sixth  Lowest  Inversion  (m) 
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